After reading an article on the subjectivity of the USDA’s rankings of American counties by scenery and climate, the author hypothesised that counties ranked lower on this scale would experience less severe impacts from the warming climate. This was based on the USDA’s criteria for the ranking, which favoured warmer counties located by the sea or near mountains–both geographic features correlated with natural disasters driven by climat change. To test this hypothesis, the author merged the dataset provided by the USDA with one gathered by the journalism nonprofit ProPublica. ProPublica, in coordination with the Rhodium Group, scored counties based on increases in the following metrics:

  1. Wet bulb temperature (mix of heat and humidity)
  2. Overall heat
  3. Loss of crop yield
  4. Sea level rise
  5. Very large fires
  6. Overall economic damage.

The scores for these metrics were added to create an overall liveability metric, and a regression performed against the counties’ subjective attractiveness.

New names:
* `1993` -> `1993...6`
* `1993` -> `1993...7`
* `January, 1941-70` -> `January, 1941-70...8`
* `January, 1941-70` -> `January, 1941-70...9`
* `July, 1941-70` -> `July, 1941-70...10`
* ...
Registered S3 methods overwritten by 'dbplyr':
  method         from
  print.tbl_lazy     
  print.tbl_sql      
── Attaching packages ───────────────────────────────────────────────────────────────────────────────── tidyverse 1.3.0 ──
✓ ggplot2 3.3.2     ✓ purrr   0.3.4
✓ tibble  3.0.4     ✓ dplyr   1.0.2
✓ tidyr   1.1.2     ✓ stringr 1.4.0
✓ readr   1.4.0     ✓ forcats 0.5.0
── Conflicts ──────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
x dplyr::filter() masks stats::filter()
x dplyr::lag()    masks stats::lag()

Call:
lm(formula = Liveability ~ Attractiveness.USDA, data = self$df)

Residuals:
     Min       1Q   Median       3Q      Max 
-14.7215  -4.1100  -0.6876   3.9188  20.5890 

Coefficients:
                    Estimate Std. Error t value Pr(>|t|)    
(Intercept)         20.81357    0.09725  214.01   <2e-16 ***
Attractiveness.USDA  0.35990    0.04269    8.43   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 5.419 on 3104 degrees of freedom
  (7 observations deleted due to missingness)
Multiple R-squared:  0.02238,   Adjusted R-squared:  0.02207 
F-statistic: 71.07 on 1 and 3104 DF,  p-value: < 2.2e-16

The hypothesis was not born out. There was a weak positive correlation, suggesting that more attractive counties were slightly more like to experience economic damage. This may be due to variance in the types of damage, so that less attractive counties had higher risk for features not correlated with attractiveness. To explore this possibility, the correlation of attractiveness with each of the variables making up overall climate risk, was calculated.

New names:
* `1993` -> `1993...6`
* `1993` -> `1993...7`
* `January, 1941-70` -> `January, 1941-70...8`
* `January, 1941-70` -> `January, 1941-70...9`
* `July, 1941-70` -> `July, 1941-70...10`
* ...
                        [,1]
Heat              0.18193631
Wet.Bulb         -0.22886514
Farm.Crop.Yields -0.08476726
Very.Large.Fires  0.50741478
Economic.Damages  0.08028562

The results somewhat bore this hypothesis out: more attractive counties were more likely to experience heat and fire-related damage, while less attractive counties were more at risk from wet bulb temperature increases (meaning higher and more dangerous humidity) and loss of crop yields. Thus, a different picture emerges than that which was originally hypothesised: less attractive counties, most of them in the center of the country, face different challenges from climate change–but not necessarily less severe ones.

To quantify these differences, counties were mapped by colour, according to the highest source of risk. As hypothesised, more and less physically attractive regions faced different risks. The Midwest and inland Northeast, both of which scored poorly on the attractiveness scale, are more likely to be affected by increases in wet bulb temperatures and overall economic losses. The Rocky Mountains and Pacific coast, rated as attractive by the USDA, faced greater risk from fire damage and sea level increases. Note that Alaska, Hawaii and Puerto Rico were not surveyed; the risk shown for these areas is a visual artifact.

Works cited

  1. USDA Natural Amenities Scale. Retrieved from https://www.ers.usda.gov/data-products/natural-amenities-scale.aspx
  2. New Climate Maps Show a Transformed United States. Retrieved from https://projects.propublica.org/climate-migration/
  3. Cartographic Boundary Files - Shapefile. Retrieved from https://www.census.gov/geographies/mapping-files/time-series/geo/carto-boundary-file.html
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQpBZnRlciByZWFkaW5nIGFuIGFydGljbGUgb24gdGhlIHN1YmplY3Rpdml0eSBvZiB0aGUgVVNEQSdzIHJhbmtpbmdzIG9mIEFtZXJpY2FuIGNvdW50aWVzIGJ5IHNjZW5lcnkgYW5kIGNsaW1hdGUsIHRoZSBhdXRob3IgaHlwb3RoZXNpc2VkIHRoYXQgY291bnRpZXMgcmFua2VkIGxvd2VyIG9uIHRoaXMgc2NhbGUgd291bGQgZXhwZXJpZW5jZSBsZXNzIHNldmVyZSBpbXBhY3RzIGZyb20gdGhlIHdhcm1pbmcgY2xpbWF0ZS4gVGhpcyB3YXMgYmFzZWQgb24gdGhlIFVTREEncyBjcml0ZXJpYSBmb3IgdGhlIHJhbmtpbmcsIHdoaWNoIGZhdm91cmVkIHdhcm1lciBjb3VudGllcyBsb2NhdGVkIGJ5IHRoZSBzZWEgb3IgbmVhciBtb3VudGFpbnMtLWJvdGggZ2VvZ3JhcGhpYyBmZWF0dXJlcyBjb3JyZWxhdGVkIHdpdGggbmF0dXJhbCBkaXNhc3RlcnMgZHJpdmVuIGJ5IGNsaW1hdCBjaGFuZ2UuClRvIHRlc3QgdGhpcyBoeXBvdGhlc2lzLCB0aGUgYXV0aG9yIG1lcmdlZCB0aGUgZGF0YXNldCBwcm92aWRlZCBieSB0aGUgVVNEQSB3aXRoIG9uZSBnYXRoZXJlZCBieSB0aGUgam91cm5hbGlzbSBub25wcm9maXQgUHJvUHVibGljYS4gUHJvUHVibGljYSwgaW4gY29vcmRpbmF0aW9uIHdpdGggdGhlIFJob2RpdW0gR3JvdXAsIHNjb3JlZCBjb3VudGllcyBiYXNlZCBvbiBpbmNyZWFzZXMgaW4gdGhlIGZvbGxvd2luZyBtZXRyaWNzOgoKMS4gV2V0IGJ1bGIgdGVtcGVyYXR1cmUgKG1peCBvZiBoZWF0IGFuZCBodW1pZGl0eSkKMi4gT3ZlcmFsbCBoZWF0CjMuIExvc3Mgb2YgY3JvcCB5aWVsZAo0LiBTZWEgbGV2ZWwgcmlzZQo1LiBWZXJ5IGxhcmdlIGZpcmVzCjYuIE92ZXJhbGwgZWNvbm9taWMgZGFtYWdlLgoKVGhlIHNjb3JlcyBmb3IgdGhlc2UgbWV0cmljcyB3ZXJlIGFkZGVkIHRvIGNyZWF0ZSBhbiBvdmVyYWxsIGxpdmVhYmlsaXR5IG1ldHJpYywgYW5kIGEgcmVncmVzc2lvbiBwZXJmb3JtZWQgYWdhaW5zdCB0aGUgY291bnRpZXMnIHN1YmplY3RpdmUgYXR0cmFjdGl2ZW5lc3MuCgpgYGB7cixlY2hvPUZBTFNFfQpSZWdyZXNzb3I8LUNsaW1hdGUuY2hhbmdlLlVTREEucmFua2luZ3M6OkxpbmVhcl9JbmZlcnJlciRuZXcoKQpSZWdyZXNzb3Ikc2ltcGxlX2xpbmVhcigpClJlZ3Jlc3NvciRsaW5lYXJfZml0cGxvdApgYGAKVGhlIGh5cG90aGVzaXMgd2FzIG5vdCBib3JuIG91dC4gVGhlcmUgd2FzIGEgd2VhayBwb3NpdGl2ZSBjb3JyZWxhdGlvbiwgc3VnZ2VzdGluZyB0aGF0IG1vcmUgYXR0cmFjdGl2ZSBjb3VudGllcyB3ZXJlIHNsaWdodGx5IG1vcmUgbGlrZSB0byBleHBlcmllbmNlIGVjb25vbWljIGRhbWFnZS4gVGhpcyBtYXkgYmUgZHVlIHRvIHZhcmlhbmNlIGluIHRoZSB0eXBlcyBvZiBkYW1hZ2UsIHNvIHRoYXQgbGVzcyBhdHRyYWN0aXZlIGNvdW50aWVzIGhhZCBoaWdoZXIgcmlzayBmb3IgZmVhdHVyZXMgbm90IGNvcnJlbGF0ZWQgd2l0aCBhdHRyYWN0aXZlbmVzcy4gVG8gZXhwbG9yZSB0aGlzIHBvc3NpYmlsaXR5LCB0aGUgY29ycmVsYXRpb24gb2YgYXR0cmFjdGl2ZW5lc3Mgd2l0aCBlYWNoIG9mIHRoZSB2YXJpYWJsZXMgbWFraW5nIHVwIG92ZXJhbGwgY2xpbWF0ZSByaXNrLCB3YXMgY2FsY3VsYXRlZC4KCmBgYHtyLCBlY2hvPUZBTFNFfQpSZWdyZXNzb3IkbGluZWFyX211bHRpdmFyaWF0ZSgpCmBgYApUaGUgcmVzdWx0cyBzb21ld2hhdCBib3JlIHRoaXMgaHlwb3RoZXNpcyBvdXQ6IG1vcmUgYXR0cmFjdGl2ZSBjb3VudGllcyB3ZXJlIG1vcmUgbGlrZWx5IHRvIGV4cGVyaWVuY2UgaGVhdCBhbmQgZmlyZS1yZWxhdGVkIGRhbWFnZSwgd2hpbGUgbGVzcyBhdHRyYWN0aXZlIGNvdW50aWVzIHdlcmUgbW9yZSBhdCByaXNrIGZyb20gd2V0IGJ1bGIgdGVtcGVyYXR1cmUgaW5jcmVhc2VzIChtZWFuaW5nIGhpZ2hlciBhbmQgbW9yZSBkYW5nZXJvdXMgaHVtaWRpdHkpIGFuZCBsb3NzIG9mIGNyb3AgeWllbGRzLiBUaHVzLCBhIGRpZmZlcmVudCBwaWN0dXJlIGVtZXJnZXMgdGhhbiB0aGF0IHdoaWNoIHdhcyBvcmlnaW5hbGx5IGh5cG90aGVzaXNlZDogbGVzcyBhdHRyYWN0aXZlIGNvdW50aWVzLCBtb3N0IG9mIHRoZW0gaW4gdGhlIGNlbnRlciBvZiB0aGUgY291bnRyeSwgZmFjZSBkaWZmZXJlbnQgY2hhbGxlbmdlcyBmcm9tIGNsaW1hdGUgY2hhbmdlLS1idXQgbm90IG5lY2Vzc2FyaWx5IGxlc3Mgc2V2ZXJlIG9uZXMuCgpUbyBxdWFudGlmeSB0aGVzZSBkaWZmZXJlbmNlcywgY291bnRpZXMgd2VyZSBtYXBwZWQgYnkgY29sb3VyLCBhY2NvcmRpbmcgdG8gdGhlIGhpZ2hlc3Qgc291cmNlIG9mIHJpc2suIEFzIGh5cG90aGVzaXNlZCwgbW9yZSBhbmQgbGVzcyBwaHlzaWNhbGx5IGF0dHJhY3RpdmUgcmVnaW9ucyBmYWNlZCBkaWZmZXJlbnQgcmlza3MuIFRoZSBNaWR3ZXN0IGFuZCBpbmxhbmQgTm9ydGhlYXN0LCBib3RoIG9mIHdoaWNoIHNjb3JlZCBwb29ybHkgb24gdGhlIGF0dHJhY3RpdmVuZXNzIHNjYWxlLCBhcmUgbW9yZSBsaWtlbHkgdG8gYmUgYWZmZWN0ZWQgYnkgaW5jcmVhc2VzIGluIHdldCBidWxiIHRlbXBlcmF0dXJlcyBhbmQgb3ZlcmFsbCBlY29ub21pYyBsb3NzZXMuIFRoZSBSb2NreSBNb3VudGFpbnMgYW5kIFBhY2lmaWMgY29hc3QsIHJhdGVkIGFzIGF0dHJhY3RpdmUgYnkgdGhlIFVTREEsIGZhY2VkIGdyZWF0ZXIgcmlzayBmcm9tIGZpcmUgZGFtYWdlIGFuZCBzZWEgbGV2ZWwgaW5jcmVhc2VzLiBOb3RlIHRoYXQgQWxhc2thLCBIYXdhaWkgYW5kIFB1ZXJ0byBSaWNvIHdlcmUgbm90IHN1cnZleWVkOyB0aGUgcmlzayBzaG93biBmb3IgdGhlc2UgYXJlYXMgaXMgYSB2aXN1YWwgYXJ0aWZhY3QuCgpgYGB7cixlY2hvPUZBTFNFfQptYXBwaW5nPC1DbGltYXRlLmNoYW5nZS5VU0RBLnJhbmtpbmdzOjpSaXNrX0lzb2xhdG9yJG5ldygpCm1hcHBpbmckbGVhZmxldF9idWlsZGVyKCkKYGBgCiMjIFdvcmtzIGNpdGVkCjEuIFVTREEgTmF0dXJhbCBBbWVuaXRpZXMgU2NhbGUuIFJldHJpZXZlZCBmcm9tIGh0dHBzOi8vd3d3LmVycy51c2RhLmdvdi9kYXRhLXByb2R1Y3RzL25hdHVyYWwtYW1lbml0aWVzLXNjYWxlLmFzcHgKMi4gTmV3IENsaW1hdGUgTWFwcyBTaG93IGEgVHJhbnNmb3JtZWQgVW5pdGVkIFN0YXRlcy4gUmV0cmlldmVkIGZyb20gaHR0cHM6Ly9wcm9qZWN0cy5wcm9wdWJsaWNhLm9yZy9jbGltYXRlLW1pZ3JhdGlvbi8KMy4gQ2FydG9ncmFwaGljIEJvdW5kYXJ5IEZpbGVzIC0gU2hhcGVmaWxlLiBSZXRyaWV2ZWQgZnJvbSBodHRwczovL3d3dy5jZW5zdXMuZ292L2dlb2dyYXBoaWVzL21hcHBpbmctZmlsZXMvdGltZS1zZXJpZXMvZ2VvL2NhcnRvLWJvdW5kYXJ5LWZpbGUuaHRtbAoK